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Abstract 

In this paper, onc-dimcnsional (ID) nonlinear wave equations 

utt - u^x + V{x)u = /(u), 

with periodic boundary conditions arc considered; V is a periodic smooth or analytic 
function and the nonlinearity / is an analytic function vanishing together with its 
derivative at u = 0. It is proved that for "most" potentials V{x), the above equation 
admits small-amplitude periodic or quasi-periodic solutions corresponding to finite 
dimensional invariant tori for an associated infinite dimensional dynamical system. 
The proof is based on an infinite dimensional KAM theorem which allows for multiple 
normal frequencies. 



1 Introduction and Results 

In the 90's, the celebrated KAM (Kolmogorov-Arnold-Moser) theory has been successfully 
extended to infinite dimensional settings so as to deal with certain classes of partial differ- 
ential equations carrying a Hamiltonian structure, including, as a typical example, wave 
equations of the form 

uu-u^. + V{x)u = f{u), f{u) = 0{u^), (1.1) 
see Wayne [^], Kuksin Q and Poschel fl^. In such papers, KAM theory for lower 



dimensional tori ||l^, |12], [0] (i.e., invariant tori of dimension lower than the number 



*The research was partly supported by NNSF of China and by the Italian CNR grant # 211.01.31. 



1 



of degrees of freedom), has been generalized in order to prove the existence of small- 
amp htude quasi-periodic solutions for (1.1) subject to Dirichlet or Neumann boundary 



conditions (on a finite interval for odd and analytic nonlinearities /). The technically 
more difficult periodic boundary condition case has been later considered by Craig and 
Wayne who established the existence of periodic solutions. The techniques used in Q 
are based not on KAM theory, but rather on a generalization of the Lyapunov-Schmidt 
procedure and on techniques by Frohlich and Spencer H. Recently, Craig and Wayne's 
approach has been significantly improved by Bourgain Q who obtained the existence 
of quasi-periodic solutions for certain kind of ID and, most notably, 2D partial differential 
equations with periodic boundary conditions. 

The technical reason why KAM theory has not been used to treat the periodic boundary 
condition case is related to the multiplicity of the spectrum of the Sturm-Liouville operator 
A = + V{x). Such multiplicity leads to some extra "small denominator" problems 

related to the so called normal frequencies. 

The purpose of this paper is to show that, allowing for more general normal forms, one can 
indeed use KAM techniques to deal also with the multiple normal frequency case arising 
in PDE's with periodic boundary conditions. 

A rough description of our results is as follows. Consider the periodic boundary problem 



for (1^) with an analytic nonlinearity / and a real analytic (or smooth enough) potential 
V. Such potential will be taken in a d-dimensional family of functions parameterized by 
a real d-vector ^, V{x) = V{x,^), satisfying general non-degenerate ("non-resonance— 
of-eigenvalue" ) conditions. Then for "most" potentials in the family (i.e. for most ^ in 
Lebesgue measure sense), there exist small-amplitude quasi-periodic solutions for ( |l . l]) 
corresponding to (i-dimensional KAM tori for the associated infinite dimensional Hamilto- 
nian system. Moreover (as usual in the KAM approach) one obtains, for the constructed 
solutions, a local normal form which provides linear stability in the case the operator A 
is positive definite. 

We believe that the technique used in this paper can be generalized so as to cover 2D wave 
equations. 

The paper is organized as follows: In section 2 we formulate a general infinite dimensional 
KAM Theorem designed to deal with multiple normal frequency cases; in section 3 we 
show how to apply the preceding KAM Theorem to the nonlinear wave equation ( |l.lD 
with periodic boundary conditions. The proof of the KAM Theorem is provided in sections 
4h-6. Some technical lemmata are proved in the Appendix. 

2 An infinite dimensional KAM Theorem 

In this section we will formulate a KAM Theorem in an infinite dimensional setting which 
can be applied to some ID partial differential equations with periodic boundary conditions. 

We start by introducing some notations. 



2 



2.1 Spaces 

For n e N, let dn G Z+ be positive even integers^. Let Z = (gineNC^": the coordinates in 
Z are given by z = {zq, zi, Z2, ■ ■ ■) with Zn = {z\, • • • , z^") E C^". Given two real numbers 
a, /3, we consider the (Banach) subspace of Z given by 

Za,p = {z e Z : \z\a,p < 00} 



where the norm \ ■ \a,p is defined as 



^ol + E l^nl^'^e"'', 



(and the norm in C^" is taken to be the 1-norm l-Zn] = X]j=i NnD- 

In what follows, we shall consider either a = and p > or a > and p = (corresponding 
respectively to the analytic case or the finitely smooth case). 

The role of complex neighborhoods in phase space of KAM theory will be played here by 
the set 

Va,p =f^X&X Za,p, 

where 1"^ is the complexification of the real torus T"^ = R'^/27rZ'^. 
For positive numbers r, s we denote by 

^a,p(r,s) = {{e,I,z) G Va,p ■■ \lra9\ < r, |/| < s'^,\z\a,p < s} (2.1) 

a complex neighborhood of T"' x {/ = 0} x {z = 0}. Finally, we denote by O a given 
compact set in with positive Lebesgue measure: ^ O will parameterize a selected 



family of potential V = V{x,S,) in (1.1) 



2.2 Functions 

We consider functions F on Da^p{r,s) x O having the following properties: (i) F is real 
for real arguments; (ii) F admits an expansion of the form 

F = J2f^z^, (2.2) 

a 

where the multi- index a runs over the set a = {ao,ai, ...) G CS'rieNN'^" with finitely many 
non-vanishing component^ a^, (hi) for each a, the function Fq, = Fa{0,I,^) is real 
analytic in the variables (9,1) G {|Im0| < r,\I\ < s^}; (iv) for each a, the dependence 
of Fa upon the parameter ^ is of class C^r{0) for some J > (to be fixed later): here 
C{y(0) denotes the class of functions which are m times differentiable on the closed set 



O in the sense of Whitney |17]. 



^We use the notations N = {0, 1, 2, • ■ ■}, Z+ = {1, 2, ■ ■ ■}. 

no no dn 

^Thus 3 no > such that = J]^ 2^" = J| J|(2^)"". 

n=0 71—0 j — 1 
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The convergence of the expansion (|2.2| ) in Da^p{r,s) x O will be guaranteed by assuming 
the finiteness of the following weighted norm: 

l|^l|z),,p(r,s),o = sup (2.3) 
where, if F« = FkiaiO^^e'^'''^^ , is defined as 

\\Fj^J2\^kia\osV\^, max |^^|, (2.4) 

k,i IpI-'^^ 

(the derivatives with respect to ^ are in the sense of Whitney). 

The set of functions F : Da^p{r, s)xO ^ C verifying (i)-^(iv) above with finite || • {r,s),o 
norm will be denoted by J^d^ p{r,s),o- 

2.3 Hamiltonian vector fields and Hamiltonian equations 

To functions F S J^£)^ p{r,s),Oj associate a Hamiltonian vector field defined as 

where denotes the standard symplectic matrix ( ^ \ • _ 

derivatives of F are defined as the derivatives term~by-term of the series ( |2.2D defining 
F. The appearence of the imaginary unit is due to notational convenience and will be 
justified later by the use of complex canonical variables. 

Correspondingly we consider the Hamiltonian equations^ 

= Fi, i = -Fe, Zn = Ud„F,„, neN. (2.5) 

A solution of such equation is intended to be just a map from an interval to the donaain 
of dehnition of F, Da^p{r, s), satisfying (^^). 

Given a real number a, we define also a weighted norm for Xp by letting^ 

WXpr''' , = (2.6) 



Notational Remark In what follows, only the indices r, s and the set O will change while 
a,a,p will be kept fixed, therefore we shall usually denote {r s) o II^Fllr,*,©, 

-Da,p(r, s) by D{r, s) and J^Da,p(r,s),o by J^r,s,o- 



^Dot stands for the time derivatives d/dt. 

*The norm || • |[j5^ p(r.s),o for scalar functions is defined in (2. J). For vector (or matrix-valued) functions 
G : Da.p{r,s) X O ^ C™, (m < oo) is similarly defined as ||G||D„,p(r,s),o = SHi ll^i ||D„,p{r,s),o (for the 
matrix-valued case the sum will run over all entries). 
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2.4 Perturbed Hamiltonians and the KAM result 

The starting point will be a family of integrable Hamiltonians of the form 

N = {Uji0,l) + I J2iM0^n,Zn), (2.7) 

nGN 

where ^ G O is a parameter, An is a (i„ x dn real symmetric matrix and (•,•) is the 
standard inner product; here the phase space T'a,p is endowed with the symplectic form 

dn/2 

n j=l 

For simplicity, we shall take, later, a;(^) = ^. 

For each ^ E O, the Hamiltonian equations of motion for N, i.e., 

dO dl dZn 

^=0., - = 0, -^ = U,„AnZn, neN, (2.8) 

admit special solutions {9, 0, 0) {O+ut, 0, 0) corresponding to an invariant torus in Va,p- 
Consider now the perturbed Hamiltonians 

H = N + P= (u;(0, 1) + IY1 iMOzn, zn) + Pie, I, z, (2.9) 

neN 

with P € J'r,s.O- 

Our goal is to prove that, for most values of parameter ^ G O (in Lebesgue measure sense), 
the Hamiltonian H = N + P still admits an invariant torus provided ||^p|| is sufHciently 
small. 

In order to obtain this kind of result we make the following assumptions on and the 
perturbation P. 

(Al) Asymptotics of eigenvalues: There exist d G N, (5 > and b> 1 such that d„ < d for 
all n, and 

= A„ ( / ) + ^n, Bn = 0{n-') (2.10) 

where A„ are real and independent of ^ while i?„ may depend on ^; furthermore, the 
behaviour of A^'s is assumed to be as follows 



Xn = n'' + o{n^), \ \ =l + o{n-^), n < m. (2.11) 
m" — n° 



(A2) Gap condition: There exists (5i > such that 

dist (cj{JdiAi), a{Jd^Aj)) > <5i > 0, Vi 7^ j ; 

denotes "spectrum of •"). 
Note that, for large i,j, the gap condition follows from the asymptotic property. 
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(A3) Smooth dependence on parameters: All entries of i?„ are cP Whitney-smooth func- 
tions of ^ with C^-norm bounded by some positive constant L. 

(A4) Non-resonance condition: 

meas{e G O : {k,u;{0){{k,u;{C)) + X{0)i{k,i^m + A(0 + KC)) = 0} = 0, (2.12) 

for each 7^ A; G Z"^ and for any A, /U G IJneN '^('^rfn^n); meas = Lebesgue measure. 

(A5) Regularity of the perturbation: The perturbation P G p{r,s),o is regular in the 
sense that ||Xp||"'^ , ^ „ < 00 with a > a. In fact, we assume that one of the following 

Oa,p(r,s),0 '-' 

holds: 

(a) p > 0, a > a = 0; (6) p = 0, a > a > 0, 

(such conditions correspond, respectively, to analytic or smooth solutions). In the case of 
d = 1 (i.e., the periodic solution case) one can allow a = a. 

Now we can state our KAM result. 



Theorem 1 Assume that N in ^2. Tj j satisfies (Al) - (A4) and P is regular in the sense 
of (A5) and let 7 > 0. There exists a positive constant e = e{d, d, b,6,6i,a — a, L, 7) such 
that if llXpll"''' < e, then the following holds true. There exists a Cantor set O-y C O 

with meas(C' \ O-y) — > as 7 — > 0, and two maps (real analytic in 6 and Whitney smooth 
in^eO) 

^■.T'^xO^^ Da,p{r, s) C Va,p, Cj:0^^ R\ 

such that for any ^ ^ and € T'^ the curve t —>■ "^{6 + u;(^)t,^) is a quasi-periodic 
solution of the Hamiltonian equations governed by H = N -\- P. Furthermore, ^(T'^,^) is 
a smoothly embedded d-dimensional H -invariant torus in Va,p- 

Remarks (i) For simplicity we shall in fact assume that all eigenvalues Aj of An are 
positive for all n's. The case of some non positive eigenvalues can be easily dealt with at 
the expense of a (even) heavier notation. 

(ii) In the above case (i.e. positive eigenvalues). Theorem 1 yields linearly stable KAM 
tori. 

(iii) The parameter 7 plays the role of the Diophantine constant for the frequency ui in 
the sense that, there is r > such that V/c G Z'^\{0}, 

^^''^^>2^- 

Notice also that is claimed to be nonempty and big only for 7 small enough. 

(iv) The regularity property o > o is needed for the measure estimates on 0\0.y. As we 
already mentioned, such regularity requirement is not necessary for periodic solution case, 
i.e., d = 1. Thus the above theorem applies immediately to the construction of periodic 
solutions for nonlinear Schrddinger equations. 

(v) . The non-degeneracy condition ( ^.121) (which is stronger than Bourgain's non-degenerate 
condition |^] but weaker than Melnikov's one ||T^) covers the multiple normal frequency 
case: this is the technical reason that allows to treat PDE's with periodic boundary con- 
ditions. 
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3 Application to ID wave equations 

In this section we show how Theorem [l| imphes the existence of quasi-periodic solutions 
for ID wave equations with periodic boundary conditions. 

Let us rewrite the wave equation (^]^) as follows 

utt + Au = f{u), Au = -u^x + y{x, Ou, x,t £R, 

u{t,x) = u{t,x + 2it), ut(t,x) = ut{t,x + 2tt), (3-1) 

where V{-,^) is a real-analytic {or smooth) periodic potential parameterized by some 
G R*^ (see below) and f{u) is a real-analytic function near u = with /(O) = /'(O) = 0. 

As it is well known, the operator A with periodic boundary conditions admits an or- 
thonormal basis of eigenfunctions (pn G L^(T),n € N, with corresponding eigenvalues /i„ 
satisfying the following asymptotics for large n 

fJ'2n-i, ^J'2n = + ^ [ V{x)dx + 0{n'^). 

For simplicity, we shall consider the case of vanishing mean value of the potential V and 
assume that all eigenvalues are positive: 

V{x)dx = , = > , V n. (3.2) 

Following Kuksin and Bourgain |^], we consider a family of real analytic (or smooth) 
potentials V{x,^), where the d-parameters ^ = (^i, • • • , ^rf) G O C R*^ are simply taken to 
he a given set of d frequencies A„. = JThTi- 



S,i = ^/a^ = A„, , i = l,--- ,d (3.3) 
where Hm are (positive) eigenvalues of^ A. 

We may also (and shall) require that there exists a positive 6i > such that 

\^J'k - i^h\ > <5i , (3.4) 

for all k > h except when k is even and h = k — 1 (in which case /i^ and /i/j might even 
coincide). 

Notice that, in particular, having d eigenvalues as independent parameters excludes the 
constant potential case V = constant (where, of course, all eigenvalues are double: ^2j-i = 
= 3^ + ^)- III fact, this case seems difficult to be handled by KAM approach even 
in the finite dimensional case. Such difficulty does not arise, instead, in the remarkable 
alternative approach developed by Craig, Wayne and Bourgain Q, Q. 

Equation (3J) may be rewritten as 

u = v, v + Au = f{u), (3.5) 
^Plenty of such potentials may be constructed with, e.g., the inverse spectral theory. 
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which, as well known, may be viewed as the (infinite dimensional) Hamiltonian equations 
ii = H^, V = —Hu associated to the Hamiltonian 

H=l{v,v) + liAu, u)+ I g{u) dx, (3.6) 

where 5' is a primitive of (— /) (with respect to the u variable) and (•, •) denotes the scalar 
product in L^. 

As in 0, we introduce coordinates q = {qq, Qi, • "),P = {Po,Pi, " " ") through the relations 

where^ A„ = y^/U^. System ( p. 51 ) is then formally equivalent to the lattice Hamiltonian 
equations 

qn = KPn, Pn = -Kqn-^, G = gij^ -^4>n)dx , (3.7) 

corresponding to the Hamiltonian function H = J2n£N ''^niQn~^Pn)^G{q). Rather than dis- 
cussing the above formal equivalence, we shall, following use the following elementary 
observation (proved in the Appendix): 

Proposition 3.1 Let V be analytic (respectively, smooth), let I be an interval and let 

te I ^ {q{t),p{t)) = (^{qn{t)}n>0, {Pn{t)}n>o] 

be an analytic (respectively, smootf^ solution of ^3. % ) such that 



sup {\qn{t)\ + \pn{t)\) n'^ e^'P < oo (3.8) 

neN 

for some p > and a = (respectively, for p = and a big enough). Then 

u{t,x) = ^ 



is an analytic (respectively, smooth) solution of ( [3. 

Before invoking Theorem 1 we still need some manipulations. We first switch to complex 
variables: Wn = ^^(^n + ''^Pn)-,Wn = :^iln - iPn)- Equations ( |37| ) read then 

dG dG 
= -iAnWn - i^^, = iA„t()„ + i-^ — , (3.9) 

OWn OWn 

where the perturbation G is given by 

G{w) = f 5(E ^^^^^Mdx. (3.10) 

^Recall that, for simplicity, we assume that all eigenvalues iJ,„ are positive. 
^Regularity refers to the components q„ and p„. 
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Next we introduce standard action-angle variables (0, /) = ((^i, • • • , Od), (Ii, • • • , Id)) in the 
{wni , • • • , ,Wn^,--- , )-space by letting, 



so that the system (3^) becomes 
iVj + Pi^ 



dt 

dWn 

dt 



dt 



dw„ 
~dt 



iXnWn + lPw^, n ^ ni,n2, - ■ ■ ,nd, (3.11) 



where P is just G with the {wm, • • • ,Wn^,Wni, • • • , variables expressed in terms of 
the (9,1) variables and the frequencies uj = (uji, ...,ujd) coincide with the parameter ^ 
introduced in (|3.3|) : 



i^i = ^,1 = 



(3.12) 



The Hamiltonian associated to ( 3.111 ) (with respect to the symplectic form dl /\ dO + 
iEn dwn A dwn) is given by 



H={UJ,I)+ ^ XnWnWn + Pi9,I,W,W,^), 
njtni,---,na 



(3.13) 



Remark Actually, in place of H in (|3.13| ) one should consider the linearization of H around 
a given point Iq and let I vary in a small ball B (of radius < s <C |/o|) in the "positive" 
quadrant {Ij > 0}. In such a way the dependence of H upon / is obviously analytic. For 
notational convenience we shall however do not report explicitly the dependence of H on 
lo. 

Finally, to put the Hamiltonian in the form (|2.9[) we couple the variables {wn,Wn) corre- 
sponding to "closer" eigenvalues. More precisely, we let Zn = {'W2n-i,W2n,W2n-i,'W2n) for 
large^n, say n > n > Tifi and denote by zo = ({'W^n} o<n<n ,{wn} o<n<n ) the remaining 



0<n<n 



conjugated variables. The Hamiltonian ( 3.13| ) takes the form 



H={u;,I) + -J2{AnZn,Zn)+P{0,I,Z,O, 



(3.14) 



where 



Diag(A2n-l, X2n, X2n-l,X2n) 
( 

/ n \ 



h 
h 



+ 



V 







A2n-1 — A2ra 





h 
h 

A2n-1 — A2n 







for n > Ud, while = Diag({An}, {A„};1 < n < n^, n / m, • • • , n^) 
dQ = n + 1 — d. 

















0^ 
















with 



The perturbation P in ( 3.14 ) has the following (nice) regularity property. 

**Compare (Al). 
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Lemma 3.1 Suppose that V is real analytic in x (respectively, belongs to the Sobolev space 
H^{T) for some k £ N). Then for small enough p > (respectively, a > 0), r > and 
s > one has 

ll^^lC':),o=o(kllp); (3.15) 

here the parameter a is taken to be (respectively, the parameter p is taken to be 0). 

A proof of this lemma is given in the Appendix. In fact, Xp is even more "regular" (a 
fact, however, not needed in what follows): ( ^.15 ) holds with 1 in place of 1/2. 




The Hamiltonian (3.14) is seen to satisfy all the assumptions of Theorem 
4, n > 1; do = n + 1 — d; d = max{do, 4}; 6 = 1; (5 = 2; 5i choosen as in (| 
Thus Theorem |l| yields the following 

Theorem 2 Consider a family of ID nonlinear wave equation ( ^g. j| j parameterized by 
^ = uj £ O as above with V{-,£,) real-analytic (respectively, smooth). Then for any < 
7^1, there is a subset of O with meas(0\C'^) ^ as 7 ^ 0, such that ( |3.1| )^p^^ 
has a family of small- amplitude (proportional to some power ofj), analytic (respectively, 
smooth) quasi-periodic solutions of the form 

u{t, x) Un{uj[t, • • • , io'^t)(t)n{x) 

n 

where u„ : T'^ ^ R and io[, - ■ ■ ,lo'^ are close to uji, • • • , uJd- 

Remark As mentioned above, our KAM theorem (which applies only to the case that 
not all the eigenvalues are multiple^ and under the hypothesis that all p„'s are positive) 
implies that the quasi-periodic solutions obtained are linearly stable. In the case that all 
the eigenvalues are double (as in the constant potential case), one should not expect linear 
stability (see the example given by Craig, Kuksin and Wayne |5|). We also notice that, 
essentially with only notational changes, the proof of the above theorem goes through in 
the case that some of the eigenvalues are negative. 



4 KAM step 

Theorem 1 will be proved by a KAM iteration which involves an infinite sequence of change 
of variables. 

At each step of the KAM scheme, we consider a Hamiltonian vector field with 



where Ni, is an "integrable normal form" and Pi, is defined in some set of the formf^ 

D{s^,r^) X O^. 

We then construct a map0 

: D{s^+i,r^+i) X O^+i c D{r^,s^) x ^ D{r^,s^) x 

^ Recall that we require that the torus frequencies are independent parameters. 
^"^Recall the notations from Section 2. 
Recall that the parameters a, p and a are fixed throughout the proof and are therefore omitted in the 
notations. 
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so that the vector field Xi{^^o<s>^ defined on I?(r,y+i, Sj^+i) satisfies 

\\Xh„o'S>^ - -'^Af^+i l|rv+l,s^+i,Cl^+l < C 

with some new normal form A/'jy+i and for some fixed i/- independent constant k > 1. 
To simplify notations, in what follows, the quantities without subscripts refer to quantities 
at the step, while the quantities with subscripts + denotes the corresponding quantities 
at the {u + 1)**^ step. Let us then consider the Hamiltonian 

H = N + P = e+{LJ,I) + ^ J2{^nZn, Zn) + P, (4.1) 

new 

defined in D{r,s) x O; teh A^s are symmetric matrices. We assume that ^ G O satisfies]^ 
(for a suitable r > to be specified later) 



|(A:,^)-i| < L^, \\{{k,u;)Id,+A,Jdr'\\ < (— )' 
7 7 

\\{{k,u;)Id,d, + iAJ,^) Id, - ® {Jd,A,)r^\\ < {\^f, (4.2) 
We also assume that 

max — - — \\ < L, (4.3) 

on O, and 

\\Xp\\r,s,o < e. (4.4) 

We now let < r+ < r, and define 

1 1 _^ 4 

s+ = -se3, e+ = 7 '^r(r — r+)e3, (4.5) 



where 



r(t) = sup^/"e-i"* 

M>1 



for t > 0. Here and later, the letter c denotes suitable (possibly different) constants that 
do not depend on the iteration stepf^. 

We now describe how to construct a set 0+ C O and a change of variables $ : x C+ = 
-D(r_|_, s+) X — > D(r,s) x O, such that the transformed Hamiltonian = + = 
H satisfies all the above iterative assumptions with new parameters s+ , e+ , r+ , 7+ , 
and with ^ G O^. 



The tensor product (or direct product) of two m x n, k x I matrices A = {aij),B is a (mfc) x (ni) 
matrix defined by 

CanB ■ ■ ■ ai„B \ 

a„i_B • ■ • amnB / 

\\ ■ II for matrix denotes the operator norm, i.e., ||Af|| = supj^^j^j^ \My\. RecaU that u and the Ai's depend 
on ^. 

^^ActuaUy, here c = J^r + (?r + + 1. 
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4.1 Solving the linearized equation 

Expand P into the Fourier- Taylor series 

k,l,a 

where k G Z*^, I E N'^ and a G (8>„gNN'^" with finite many non-vanishing components. 
Let R be the truncation of P given by 

k,\l\<l 

+ E ^'wae'<'='^> + E PrnJ^^'^'^z^, (4.6) 

fc,|a|=l fc,|«|=2 

with 

2\l\ + \a\ = 2 E ^j+El"il-^- 

It is convenient to rewrite R as fohows 

R(ej,z) = E ^fe/oe'^'-'^ 
fc,|i|<i 

+ E(^''^^)e'^'''^ +E(4^-^^>'^'''^ (4-7) 

k,i k,i,j 

where RfjRj^ are respectively the di vector and {dj x di) matrix defined by 

. I g.-.<M„,|^^„^„. . i±l I WL,-,<«),,|„„,.„, ,4,8) 

Note that i?,'^^- = (R'^if . 

Rewrite H as H = N + R + (P — R). By the choice of s+ in ( ^.51 ) and by the definition 
of the norms, it follows immediately that 

\\XR\\r,s,0 < \\Xp\\r,s,0 < £• (4.9) 

Moreover s+, e+ are such that, in a smaller domain D{r, s+), we have 

\\Xp_R\\r,s+ < ce+. (4.10) 

Then we look for a special F, defined in domain = D{rj^, s+), such that the time one 
map (l)^p of the Hamiltonian vector field Xp defines a map from Dj^ — > D and transforms 
H into 

More precisely, by second order Taylor formula, we have 
Ho<i?p = {N + R) o + {P - R) o <i)\ 
= N + {N,F]+R 

+\ I ds ['{{N + R, F}, F} o (t)'pdt + {R, F} + {P-R)o c/,),. 
2 JO Jo 

= N+ + P+ 

+ {N,F} + R-PooO-{uj'j)-Y.iPnn^n,Zn), (4.11) 

neN 
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where 



^' - j -Qjde\i=Q^z=o, Rnn - j 



17=0,^=0, 



P+ = \ C ds ['{{N + R, F}, F} o X'pdt + {R, F} + {P - R) o 0],. 
^ Jo JO 



"1 rs 

ds 
10 Jo 

We shall find a function F of the form 



F{e,I,z) = Fo + F,+F2= i^Hoe^<'^''^/' + E(^''^*)e'^'''^ 

|«|<i,|fc|^o «eN 

+ {F^z.,z,)S''''K 
|fc|+|i-jVo 

satisfying the equation 

{iV, F}+R- PoOO - {OJ', i^nn^n, Zn) = 0. 



(4.12) 



(4.13) 



Lemma 4.1 Equation ( j.l^ is equivalent to 



Fkio = ii{k,Lo))-^PMo, A;y^O,|/|<l, 
i{k,u;)I,^ + A,MF^^-Ft^mA,) = iRl 



+ |i-i|^0. (4.14) 



Proof. Inserting F, defined in ( [4.12 ), into ( [4.13| ) one sees that ( [4.13 ) is equivalent to the 
following equation^ 



{7V,Fo} + Po-(w',/) =0, 
{Ar,Fi} + Pi = 0, 

{N, F2} +P2-Y. i^nnZn, Zn) = 0. 



(4.15) 



neZ 



The first equation in ( 4.15| ) is obviously equivalent, by comparing the coefficients, to the 
first equation in ([4.14[) . To solve {N, Fi} + Pi = 0, we note that^ 

{N,Fi} = {diN,deFi) + {V,N,JV,Fi) 

= {diN,deFi)+Y{'^z,N,UaV,^Fi) 



iYii {k,u;)Ft,z,) + {Az„J,^Ft))e'^'^''^ 

k,i 

iYm,^)Id,+AiWFtzi)e'^'^^'^ . 



(4.16) 



kA 



Recall the definition of Pi in ([4.6[) 
Recall the definition of TV in (|4.lD 



13 



It follows that Fj' are determined by the linear algebraic system 
i((A;, u;)Id^ + AJdjFt + = 0, i G N, A; G 

Similarly, from 
{N,F2} = {diN,deF2) + Y.{V,^N,iJdy,,F2) 

i 

= i {{{k,u;)F^,Zi,z,) + {AiZ,,Jd^{F^^^Zj) + {A^zj,Jd^F^,z,))e'(^^^^ 

= i {{{kMF^^^^.z,) + {{A,Ja^Ff,-Fp,^A,)z,,z,))e^^^'') 

|fc| + |i-j|^0 

= i m,u;)F^,+A,J,^F^,-F^,J,^A,)z„z,)e'^'''^^ 



(4.17) 



it follows that, Fj- is determined by the following matrix equation 

{{k,u;)^ + A,^)F^, - F^,{JdM = i^i, 1^1 + |i " j| / 



(4.18) 



where F^^, Rj^ are dj x di matrices defined in ( 4.12| ) and ([4.7|). Exchanging i,j we get the 
third equation in ( [4.14 ). ■ 

The first two equations in ( 4.14| ) are immediately solved in view of (^]^). In order to solve 
the third equation in (4.14), we need the following elementary algebraic result from matrix 
theory. 

Lemma 4.2 Let A,B,C be respectively n x n,m x m,n x m matrices, and let X be an 
n X m unknown matrix. The matrix equation 

AX- XB = C, (4.19) 

is solvable if and only if 1^ ® A — B ® In is nonsingular. Moreover, 

\\X\\ < \\{I„,®A-B^In)~^\\ ■ \\C\\. 

In fact, the matrix equation ( 4.19| ) is equivalent to the (bigger) vector equation given by 
{I ® A — B ® I)X' = C where X', C are vectors whose elements are just the list (row by 
row) of the entries of X and C. For a detailed proof we refer the reader to the Appendix 
in [||] or 0], p. 256. 

Remark. Taking the transpose of the third equation in ( 4.14| ), one sees that [F^^'^ 

satisfies the same equation of F^^. Then (by the uniqueness of the solution) it follows that 

pk ^ (pk\T 
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4.2 Estimates on the coordinate transformation 

We proceed to estimate Xp and (^\p. We start with the following 
Lemma 4.3 Let Di = D{ls, r+ + {{r - r+)), < i < 4. Then 

\\Xf\\d,,o < cj-'T{r-r+)e. (4.20) 



Proof. By ( |4.2D , Lemma LI and Lemmata 7A, 7^ in the Appendix, we have 



WFtWo = mk,u;)h^+AM~'R^\\ < \\{{k,u;)h,+AJ,r'\\ ■ \\Ri\\ 

WFt^Wo < \\i{k,io)h,d,+{A.M^Id,-Id,^iJd,A,)r'\\-\\R'^^\\ 

< c-f-mR'ljl \k\ + \i-j\^0. (4.21) 

Where || • Ho for matrix is similar to (12.41) . 
It follows that 

^WFoWo < ^(E l/«o| • |/'| • \k\ ■ |e'<^'^> I + E l^'l • h\ ■ \k\ ■ W^'''^ I 

< c j-'T{r-r+)\\Xn\\ 

< c j-''T{r-r+)e, (4.22) 
where T{r — r+) = sup;. \k\'^e~^^^i^^~'^'^\ Similarly, 

II^/||d.,o = E \Rkio\ ■ |e'<'^'^> I < c7-^r(r-r+)e. 

Now we estimate ||Xpi H/jj^^. Note that 



e 



\D2,0 



<c ^-TEl^'lel'^l'^ < C7'T||-^||. (4.23) 



It follows that 



< c J-'^Ty \\—!-\\i''e'P < C7-Te, 



by the definition of the weighted norm. 
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Note that 



I f2 II 



< c 7~=r 



E(4 + (Fl^ f)^^^^'''^ Wo 



dzi 



Similarly, we have 

\\Xp2\\D2,o < c7"'Te. 
The conclusion of the lemma follows from the above estimates. 



(4.24) 
(4.25) 



In the next lemma, we give some estimates for (p^p. The following formula ([4.26|) will be 



used to prove that our coordinate transformations is well defined. (4.27) is for proving the 
convergence of the iteration. 



Lemma 4.4 Let r] = es ,Dj^ = D{r^ + — r+)A7]s),i = 1,2. We then have 



if e <^ {^^ ^)2. Moreover, 



\\D<j)'p-Id\\D, < c7-=re. 

3'' 



(4.26) 



(4.27) 



Proof. Let 



\D"'F\\d,o = max{| 



Q\i\ + \l\+P 

WdFd^ 



F\D,oAi\ + \l\ + |a| = m > 2}. 



Note that F is polynomial in / of order 1, in z of order 2. FromP^ (|4.2E| ) and the Cauchy 
inequality, it follows that 

\\D'^F\\d,,o < cj-'Te, (4.28) 

for any m > 2. 

To get the estimates for (pp, we start from the integral equation. 



id+ j Xp o (j)p ds 
Jo 



4>F ■ -^1,7 Dr,, < t < 1, follows directly from (lOSl) . Since 

D(j)^p = Id+ [^{DXp)D(l)'pds = Id+ J{D^F)D(j)'pds, 
Jo Jo 

it follows that 

\\D(l)^p - Id\\ <2\\D^F\\ < C7-Te. 
The estimates of second order derivative D'^(j)p follows from (4.28). 



(4.29) 



^® Recall the definition of the weighted norm in (2 
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4.3 Estimates for the new normal form 



The map (jip defined above transforms H into H+ = N+ + P+(see (pl|) and (|4.13D ) with 



N+ = e+ + {uj+,y) + \Y. ^^t^u z^) (4.30) 



where 

e+ = e + Pooo, uJ+ = uJ + Pm{\l\ = 1), ^+ = + 2RI. (4.31) 

Now we prove that shares the same properties with N . By the regularity of Xp and 
by Cauchy estimates, we have 

\uj+-uj\<e, \\R%\\<ei~^ (4.32) 
with (5 = a — a>0. It follows that 

provided I A; I '^^e < c{'^^-^j^). 
Similarly, we have 

||((A;,.. + Pm^)Id,d, + {AtJd,) ® Id, - Id, ® {Jd,Ap)-^\\ < {\^f, (4.34) 



provided |A;|'^'^e < c (7*^ — 7^). This means that in the next KAM step, small denomi- 
nator conditions are automatically satisfied for \k\ < K where K'^'^e < c (7"^ — 74!^). 
The following bounds wil be used later for the measure estimates. 



I Q^i \o<^^ l^-^^r^lc^ < (4-35) 
for |Z| < (? by the definition of the norm. 

4.4 Estimates for the new perturbation 

To complete the KAM step we have to estimate the new error term. 
By the definition of (\^p and Lemma 4.4, 



is well defined in -Di^. Moreover, we have the following estimates 



< [[-'^{{Af+iJ.Fl.FlIlD,, + ll-'^P-ijIlD,, 

< c 7-=r2et < ce+ (4.36) 
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by (|4.90 and Lemma fL. 
Thus, there exists a big constant c, independent of iteration steps, such that 



\\Xp^\\r+,s+ 

The KAM step is now completed. 



2" 



(4.37) 



5 Iteration Lemma and Convergence 

For any given s,e, r, 7, we define, for al\ v > 1, the following sequences 



4 

7(1 -E 2-^) 

i=2 

1 i ^ _ ^ 

1 ''""^ 



7t/ 



i=0 



(^;ii(7fi 



(5.1) 



where c is the constant in ( 4.37 ). The parameters ro, eo, 70, -^>0) -sO) Kq are defined respec- 
tively to be r, e, 7, L, s, 1. 



Note that 



^W=ri[r(r.-i-r,)]2(|)^ 



is a well defined finite function of r. 



5.1 Iteration Lemma 

The analysis of the preceeeding section can be summarized as follows. 

Lemma 5.1 Suppose that eo = e{d,d,5,5i,a — a,L,T,'y) is small enough. Then the fol- 
lowing holds for all u >0. Let 
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be a normal form with parameters ^ satisfying 



\{k,u;,)-'\ < \\{i{k,LO,)Ia,+A';jaJ-'\\ < (^) 

7i/ lu 



mk,u;.)Id.,^ + (A^J,J 7,^- - ® (Jrf,A,^))-1 < (^)'^" (5.2) 

II' 

on a closed set Ou of for all k ^ 0,i, j E Z. Moreover, suppose that uji^{^) , {^) are 
smooth and satisfy 

on Oi^ in Whitney's sense. 
Finally, assume that 

Then, there is a subset Oi,+i C Oi,, 



where 



a+l = a-U|fc|>x.+l7^^+■'(7.) 



|(fc,c^„+i>-i|>^, ||((fc,^,>72„+(A^-+ij,j-i||>(M_)d,or _^ 
C G C?^ : I ||((fc,c^.+i>/d.d.+(^^-+iJd.)0/d.-/d.0(Jd.^^-+i))-i||>(^)'='' 



wii/t Wjy+i = LOn + Pqiq, and a symplectic change of variables 

: X ^ D„ (5.3) 

such that H^^i = o defined on D^+i x Ojy+i, /las i/ie form 

H^+i = e^+i + (a;^+i, /) + ^(AJ'+^z^, 2^) + P^+i, (5.4) 

satisfying 

, ^'(^"+'(0-4^ 1 ^ -5 

5.2 Convergence 

Suppose that the assumptions of Theorem 1 are satisfied. To apply the iteration lemma 
with u = 0, recall that 

eo = e, 70 = 7, sq = s,Lq = L, Nq = N,Po = P, 
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|{fc,a;>-i|<Hl,||((fc,^)/,^+A,j,j-i||<(Ml)5 or 
Co = U G O : I ||((fc,^>/d,d,+{A,Jd,)®/d^.-/d,®(Jd.A,))-i||<(^)^^ 



(with e and 7 small enough). Inductively, we obtain the following sequences 

= $1 o . . . o : Li^+i X ^ Do, 1^ > 0, 



Let = n'^^QOiy. As in [|15|, thanks to Lemma we may conclude that N^,'^'^, D"^" , uj^+i 
converge uniformly on Doo x Oj = D(^r, 0,0) x with 

iVoo = 600 + i^^oo,!) + {AocZ,z) = Coo + (t^oo,^) + ^{Afzi,Zi), 

Since 

e^+i = C7-T(r,, - r,,+i)e^ < (c7-"^(r)e)(t)'. 

It follows that Eiz+i ^ provided e is sufficiently small. 

Let be the flow of Xh. Since H o vj/'' = H^j^i, we have that 



The convergence of ^D"^^ ,bJy^i,XH^+i^ implies that one can take limit in ( |5.7] ) so as to 
get 

</,*,o*- = *-o0*,^, (5.8) 

on D(ir,0,0) X O^, with 

^'"^ :I?(^r, 0,0) X ^ X M'^. 

From (U) it follows that 

</,*,(^-(T'^ X {e})) = ^°°0*;v^(t" X u}) = M/-(T'^ X {e}), 

for ^ G O^. This means that ^'°°(T'^ X {w}) is an embedded torus invariant for the 
original perturbed Hamiltonian system at G O^. We remark here the frequencies u;oo(0 
associated to ^""(T'^ x {^}) is slightly different from ^. The normal behaviour of the 
invariant torus is governed by the matrix = X^vGN^f- ■ 
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6 Measure Estimates 



At each KAM step, we have to exclude the following resonant set of ^'s: 

\k\>K^,i,j 

the sets T^kjT^ki^'^kij being respectively 

and {cj G a : \\Mo^\\ > {}-^f], (6.1) 

where 

M2 ={Ku:,)h^d,+{Apd,)®h^-h,®{JdA"i)- (6-2) 

We include in the set € O : ||Af(ci;)^^|| > C} also the ^'s for which M is not invertible. 
Remind that u:,{i) = ^ + E -=0 ^000(6 withQ | E PiooiO^^ < e, A- = A + 2j:. 

with ||E.«°'1l=o(er^). 

Lemma 6.1 There is a constant Kq such that, for any i,j, and \k\ > Kq, 



Proof. As it is well known 



meas ijVQ < 777- • 



The set TZ^^ is empty if i > const |fc|, while, if i < const from Lemmata |7.6| , |7.7| there 
follows that 

meas(7^^J < c j^^j^- 

We now give a detailed proof for the most complicated estimate, i.e., the estimate on the 
measure of the set T^^jj- Rewrite A^2 as 

with 

Aij = {k, LOu+i)Id,dj + Diag(/rf^/2, -Id,/2) Id, - Kh, ® Diag (-4^/2, 4,/2)- (6.3) 

The matrix Aij is diagonal with entries X^ij = {k,(jjy) ^Xi^Xj in the diagonal where Aj, Xj 
are given in ( p.lOj ) and ± sign depends on the position. B^^ is a matrix of size 0{i^^ 
since A'^ = Ai + Bi + 0{i-^) = Ai + 0{i-^) by ( pl| ) and (|432l). 



^Recall (|^, 
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In the rest of the proof we drop in the notation the indices i,j since they are fixed. Now 
either all X^ij < \k\ or there are some diagonal elements Xkij > \k\. We first consider 
the latter case. By permuting rows and columns, we can find two non-singular matrices 
Qi, Q2 with elements 1 or such that 

where Aii,A22 are diagonal matrices and An contains all diagonal elements X^ij which 
are bigger than Moreover, defining Q3,Qi,D as 



r, _i I Mn-M -(^ii+^ii)~'^i2 



and 

D = A22 + B22 - B2M11 + Bii)-^Bi2 = A22 + 0{r^+j-^), (6.5) 

we have 



Q3Q1 [A + B-^') =y^^' ^^-'^ ^ J (6.6) 

For ^ G O such that D is invertible, we have 

Since the norm of Qi,Q2,Q3,Q4:,iAu + Bu)^^ are uniformly bounded, it follows from 
that 

{UO,:\\{A + Bn-^\\> {^-^f} C{UO^:\\D-^\\ >c(^)'^"'}. (6.8) 

If all Xkij < c I A;| we simply take D = A + B^ . Since all elements in D are of size 0(|A;|), 
by Lemma "LQ in the Appendix, we have 

{ieO,:\\D~^\\ >c (}^f]c{ieO,:\<ietD\ < c{-^f]. (6.9) 

Let N denote the dimension of D (which is not bigger than d^). Since D = A22 + 0{i-^ + 
the N^^ order derivative of det-D with respective to some is bounded away from 
zero by ^\k\^ (provided is bigger enough). Prom (|6.8|), ( |6.9[ ) and Lemma 7/?, it follows 
that 

meas7^L■ = meas^ G C : \\(,A + B'')-^]] > {\^ f} 
< meas{^ eO^:\detD\ < c {—1^)^} 

This proofs the lemma. ■ 
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Lemma 6.2 If i > c then TZ'^- = 0; //max{i,j} > c \k\'>-''-,i j for b > I or 
\i — j\ > const \k\ for b = 1, then TZ^^j = where the constant c depends on the diameter 
ofO. 

Proof. As above, we only consider the most complicated case, i.e., the case of TV^^ij- Notice 

1 

that max{i, j} > const \k\ for 6 > 1 or |i — j| > const \k\ for 6=1 implies 

|A,±A,| = (/-i^)(i + 0(r^+r^)) 

> \\j - + + 0{i-'+j-')) > const |A:|. (6.11) 

It follows that Aij defined in (|6.3|) is invertible and 

ii(A,r'ii<i^r'- 

By Neumann series, we have \\{Aij < 2|/i;|^^ for large k (say \k\ > Kq), i.e, 

= 0- ■ 

Lemma 6.3 For b > 1, we have 

meas(|J = meas |J {Til U TZli U TZlij) < C71TJ. 

iy>0 u,\k\>K„,i,j 



Proof. The measure estimates for IZ^ comes from our assumption ( p. 12 ). We then consider 
the estimate 

meas(U (J \Jni,^) , 

" \k\>K^ i,j 

which is the most complicate one. 

Let us consider separately the case b > 1 and the case 6 = 1. We first consider 6 > 1. By 
Lemmata 6J, 6^, if |A;| > Kq and i 7^ j, we have 



meas 



(U TZkij) = meas ( (J 7^f,) < c > c 
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iftj;i,j<C\k\b- 



(6.12) 



For i = j. As in Lemma p.l| , we can find Qi,Q2 so that (|6^) holds with the diagonal 
elements of An being < k^ujy > ±2Aj and A22 =< k,uj^ > I. Repeating the arguments 



in Lemma 6.1, we get ( |6.9| ) and 



C {Mdetl?| < ci-l^f} 



{C:ll\{k,u.)+Oiz-')\ < ci-^f} 



C {C:\{k,u;,)\ < c(^ + ^)} = Qt 



(6.13) 
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Since Q^j C . for i > io, using ( 6.10|) , we find that 



«07 , l_^ 

5' 



IT-I, 1 



for any io- Following Poschel ([|15|), we choose iq = {'—^ — )^+'^, so that 



meas 



{[jTZkul < c(-^)iTi (6.14) 



|A:| 



Let r > max{d + 2 + (d + 1)^ + 1}. As in (|]T|), ( |6lD , we find 



meas( (J U7^^,,(7.)) = meas ( IJ IJ 



i+s , 



+meas (U|fc|>;^-^ |J 7^^ij(7^+l)) < c K^^-f 

i 

The quantity meas(lj,y D\k\>Ki, '^kij) ^'^^^ bounded by 

Emeas( IJ \Jni,^{^,)) < c ^t^s R-' < C7^, (6.15) 

i^>l \k\>K^ i,j u>0 

provided r > max{(i + 2 + ^^j, (d + 1)^ + 1}. This concludes the proof for 6 > 1. 

Consider now 6 = 1. Without loss of generality, we assume j > i and j = i + m. Note 
that Lemma implies TZ\j = for m > C\k\. Following the scheme of the above proof, 
we find 

k,i,j k,i,m k,m<C\k\ i 

,io+m ,io+m 

). (6.16) 

k,m.<C\k\ i<io 

where 

7 1 

Qkio,io+m = : I {k, i^u) + m| < c ( j^^pT + — 5)}- 
Again, taking Zq^^ = -^^L — ^ -^^g have, for fixed A;, 

ij m<C|fc| ' ' 



< c lA^Kp^)^ . (6.17) 



As in the case 6 > 1, we have that meas(ljj^ lJ|;j,|->;^^ Uij ^fcij) is bounded by 0(71+*) if 
T>(d+ 1)^ + 1. ■ 



Remark In ( 6.13 ), IdetZ?! = K^;"^) + 0(i~'')| (guaranteed by the regularity property) 
is crucial for the proof. But it is not necessary for the periodic solution case, i.e., d = \. 
since Tlk,i+mi = if, ^ > c > |A:| ^ 1 are sufficiently large. 
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7 Appendix 

Proof of Proposition \3.]\ From the hypotheses there follows that the eigenfuctions (f>n are 
analytic (respectively, smooth) and bounded with, in particular, 

sup(|(/)^| + < const Hn ■ 

R 

Thus, the sum defining u{t,x) is uniformly convergent in / x [0,27r]. Since 

"0 ' fnT.^M<P. 



k 

one has 



knl < const , l^nl < const \n < const 



-np 



(7.1) 



e 

Or, < const -r . 

Thus (if a is big enough, in the smooth case) u(t, x) is a function and 

Utt + Au = ^ + -^Mn 

V A y An 

= E(/ fW)<l^n)K = f{u) , 

where in the last equality we used the fact that f{u) is a smooth periodic function. ■ 
Lemma 7.1 

II^G'llD(r,s) < l|-^llD{r,s)l|G'||£)(r,s)- 

Proof. Since {FG)kip = J2i Fk-k',i-i',p-p'Gk'i'p'- we have that 
l|i^G'b(,,,) = supY,\iFG)kip\ \y\' |^°|el'=l^ 

,p-p'Gk'i'p'\ |y|'|z"|e 

^ kip I' 

= l|-^llD(r,s) 11^1115(^,5) (7.2) 



and the proof is finished. 
Lemma 7.2 (Cauchy inequalities) 

\\P'9i\\D{r-a,s) < Ca^^ \\F\\ D{r,s) , 

and 
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Let {•, •} is Poisson bracket of smooth functions 



^.dFdG dFdo. dF .dG. 



where J^- are standard symplectic matrix in M*. 
Lemma 7.3 // 



then 

Proof. Note that 
_d 

Since 



\\XF\\r,s < e\\\XG\\r^s < e\ 



W II ^ — 1 — 2 / // 

I ^{i^,G} 1 1 r-<T,r;s <CCJ £6, 



+ E((^-.-n,-/d.G,J + (F,^, Jrf,G,,,J) (7.4) 



(-^6»2„, < c a ^\\F;,^\\\ ■ \\Gy\ 

.-2| 



IK^e,G,,J||^(,_,_i,) < c s-'\\Fe\\-\\G,J 
\mz^,Ge)\\r>i^r,\s) < c s-^F,J-\\Ge\\ 
\m,Gg,„)\\nir-.,s) < c ^-^F.ll • ||G,J| 

||(i^...„,Jd.G,,>|lz5(.,is) <c 5-'||i^.J|-||G.J|iV^ (7.5) 
it follows from the definition of the weighted norm(see (|2.6| )), that 



\\X{F,G}\\r-a,rjs < C a ^T] ^e'e". 

1 4 

In particular, if 77 ~ es , e', e" ~ e, we have ||r_CT,r;s ~ . 



Lemma 7.4 Let O &e a compact set in M*^ /or which (4-^) holds. Suppose that /(^) and 
uj{(,) are C" Whitney-smooth function in £ O with G^y norm bounded by L. Then 

_ /«) 

is C"* Whitney-smooth in O witl^^ 

\\g\\o < c-f-^lkl^L, 

Proof. The proof follows directly from the definition of the Whitney's differentiability. ■ 
Recall the definition in (12 
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A Similar lemma for matrices holds: 



Lemma 7.5 Let O be a compact set in R"^ for which (4-^) holds. Suppose that B{^),Ai{^) 
are C" Whitney-smooth matrices anduj{^) is a Withney- smooth function in(^ £ O bounded 
by L. Then 

C{0 = BM-\ 

is C™" Whitney- smooth with 

\\F\\o < cj-'\k\^L, 

where M stands for either {k, a;)/^. +AiJ(i- if B is [di x di)-matrix, or {k, ijj)Iii^dj + {-^iJdi)® 
Idj - Idi ® ( Jdj^j) if B is {didj X didj)-matrix, 

For a N X N matrix M = {aij), we denote by |M| its determinant. Consider M as a linear 
operator on {R^ , | • |) where = X) Let ||M|| be its operator norm. It is known ||M|| 
is equivalent to norm max | aij \ . Since a constant depends only on the space dimension 
and two fixed norms is irrelevant, we will simply denote ||M|| = max|ajj|. 

Lemma 7.6 Let M be a N x N non-singular matrix with \\M\\ < c\k\, then 
{lo : \\M-^\\ > /i} C {w : |detM| < c ' } 



Proof. Firstly, we note that if M is a nonsingular N x N matrix with elements bounded 
by \mij\ < m, by Cramer rule, the inverse of M is = -pgjadjM. Thus 



M-^ll < c 



N-l 



|detM| 

where the constant depends on N. In particular, if m = const|fc|, |DetM| > ^-^^ — then 

\\M-^\\ < ch. 

This proofs the lemma. ■ 

In order to estimate the measure of TZ'^'^^, we need the following lemma, which has 
been proven in [|l8| pH]. A similar estimate is also used by Bourgain p]. 



Lemma 7.7 Suppose that g{u) is a function on the closure I, where I C is a finite 
interval. Let Lh = {u\ \giu)\ < h}, h > 0. If for some constant d > 0, \g^"^\u)\ > d for 
all u £ I, then meas (Ih) < ch^ where c = 2(2 + 3 + • • • + m + d~^). 



For the proof of Lemma 3.1, we need the following 
Lemma 7.8 



if p < r,q Zp where C depends on r — p. 
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Lemma 7.9 

5](l + |n-i|)-^|jr < c|nr, ^ \qj\{l + \n-j\)-^\nr<C\q\a 
if K > a + 1, q £ 2a^p=o where C depends on K — a — 1. 
The proofs of the above two lemmata are elementary and we omit them. 



Proof of Lemma 3. 1 : Here we give a direct proof. It is clearly enough to consider the 
case of f{u) being a monomial u^~^^ for some > 1. From (|3lo|) , one can see that 
the regularity of G implies the regularity of G. In the following, we shall give the proof 
for G. Suppose that the potential V{x) is analytic in |Im3;| < r (respectively, belongs to 

Sobolev space H^) then the eigenfunctions are analytic in |Imx| < r (respectively, belong 
to If we let (t)i{x) = ^afe^<"'^> then(see, e.g.,[|) 

|a"| < ce"!*""!'' respectively |a"| < c (1 + |n - 

Recall that 



_1io ' ' ' QiN 



where 



N 

™o+"iH hnjv=0 s=0 

with I a"^" I < c e~l*^~"'' I'' (respectively, I a"^'' I < c {I + \ns - is\~^~'^). 
In what follows, we assume either a = 0, p > or a > 0, p = 0. Since 

Q.ii ' ' ' QiN 



it follows that 



< c 



\\G,t^,, = \\G,,\\+Y.\G,,\\jr'^^'' 

N 



.■■■■»jv. s=l 

iqH ^-n^—O 



N 



< c E (i + u-"oi)-^ijrei^'i''-i"°-^'iM n(i + 



N\^\aJj\p-\no-j\r { Y\ 4- In - i -2 ^-\ns-is\r 



npH hn^— 



< c 



< c 



E l^orel'^^l" [n(l + K - is\)-''-'e-^^'-''^'\qj) 
n, \s=l / 

qH 

^ N / ^ 

E (I E ns\re^^s^i^^\P J|(i + In, - z,!)"^" Vl^-^^lH^i, 
i.--,sjv; s=l \s=l 



"l''"'"JV 
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< c 



< c 



< c 
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